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A phase diagram, defined by the amplitude square and phase of scattering coefficients for ab¬ 
sorption cross-section in each individual channel, is introduced as a universal map on the electro¬ 
magnetic properties for passive scatterers. General physical bounds are naturally revealed based on 
the intrinsic power conservation in a passive scattering system, entailing power competitions among 
scattering, absorption, and extinction. Exotic scattering and absorption phenomena, from resonant 
scattering, invisible cloaking, coherent perfect absorber, and subwavelength superscattering can all 
be illustrated in this phase diagram. With electrically small core-shell scatterers as an example, 
we demonstrate a systematic method to design field-controllable structures based on the allowed 
trajectories in the phase diagram. The proposed phase diagram not only provides a simple tool to 
design optical devices but also promotes a deep understanding on Mie’s scattering theory. 

PACS numbers: 42.25.Bs, 81.05.Zx, 42.25.Fx, 78.67.Bf 


Even though studies on scattering of light by small 
particles can been tracked back to Mie’s seminal work 
one century ago [T] , it was realized in the past decade to 
have an efficient way in manipulating and designing elec¬ 
tromagnetic properties of light with nanostructures [2]. 
In particular, unusual scattering characteristics, such as 
resonant scattering nn, coherent perfect absorption 
HO], invisible cloaking mHis], subwavelength superscat¬ 
tering [T7HIS], and minimum-scattering superabsorbers 
[20l|2T], are revealed in a single, isotropic, and homoge¬ 
neous scatterer by embedding multi-layered structures. 
Combined with state-of-the-art fabrication technologies, 
these meta-structures sever as functional nano-devices 
with promising applications in light harvesting [22II24] . 
heat generation by metal nanoparticles [23 EH], opti¬ 
cal nanocircuits and nonlinear optical processes 

mEg. 

To have exotic scattering properties, a variety of spe¬ 
cific conditions are asked to be satisfied for the desired 
cross-sections constituted by scattering coefficients. Un¬ 
doubtedly, a better understanding in the scattering coef¬ 
ficients provides an access to design nanostructures, such 
as information about the limitation, energy assignment, 
and robustness on the corresponding extrinsic field re¬ 
sponse of real scatterers. For the working frequency 
of interests, for example, in the visible spectra many 
metals possess negative permittivities associated with a 
strong dispersion, which introduces real loss effect and 
suppresses the desired functions m- As possible mis¬ 
matching in physical parameters may happen, it is nat¬ 
urally to seek optimized invisible cloaks or performance 
boundary in a cloaked sensor with the consideration of 
intrinsic loss in reality m- 

In this Letter, we study the general relation between 
the amplitude and phase in the scattering coefficients for 
a passive spherical scatterer under any frequency exci¬ 
tation. By decomposing absorption powers for each in¬ 
dividual channel, a phase diagram is introduced by re¬ 
casting absorption cross-section in terms of the ampli¬ 


tude square and phase of scattering coefficients, which 
acts as a universal map to design passive scatterers. Not 
only all the physically allowed regions can be defined to 
satisfy the intrinsic power conservation, but also all the 
exotic electromagnetic properties in the literature can be 
illustrated in this phase diagram. Moreover, as a demon¬ 
stration, we take electrically small core-shell scatterers 
as an example to illustrate a systematic way in design¬ 
ing the composition of subwavelength-structures with re¬ 
quired scattering and absorption properties. Our work 
not only promotes an alternative interpretation of Mie’s 
scattering theory, but also provides a compact solution to 
realize a variety of functionally subwavelength-devices. 

We consider a linearly polarized plane wave with time 
evolution at the frequency cj, which is illuminat¬ 

ing on a single spherical object. The object could be 
made by multiple layers of uniform and isotropic media 
with complex permittivity and permeability. The real 
and imaginary parts of permittivity (permeability) are 
denoted as e = e'Pie" (/i = /i' + i/i"), where e" {fi") is as¬ 
sumed to be a positive real number for a passive medium. 
Without loss of generality, the surrounding environment 
is taken as non-absorptive, non-magnetic, and free of ex¬ 
ternal sources or currents, i.e., eo = 1 and fio = 1. Due to 
the free divergence in the governing Maxwell’s equation, 
the electric field E and magnetic field H in the environ¬ 
ment can be expressed by two auxiliary vector potentials, 
i.e., the transverse magnetic (TM) and transverse elec¬ 
tric (TE) modes, which are respectively generated by two 
scalar spherical wave equations [3 EH EH). Each scalar 
function can be built by an infinite series with unknown 
coefficients determined through the boundary conditions. 
By following the conventional notations, let the scatter¬ 
ing coefficients be S™ and for TM and TE modes 
in each spherical harmonic channel labeled by the in¬ 
dex n, respectively. It is known that the corresponding 
absorption and scattering cross-sections, and 
defined as the total power absorbed and scattered by a 
single scatterer with respect to the unit intensity of inci- 
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dent plane wave, can be expressed as 

oo 

^abs = + af = (1) 

n=l 

- E +sD+is^p+isr p), 

n=l 

^■,.,"g (2"+^l)A" (|gTM|. + |gTE|2)^ (2) 

n=l 

where A is the wavelength of incident wave in vacuum 
[3^ [33] . Due to the conservation of power, the ab¬ 
sorption cross-section is required to be equal or larger 
than zero for each channel, i.e., = —|^(2n + 

l)(Re{5'i™’^^^} + > 0. Moreover, for a 

given radius of particle, denoted as a, the value of size 
parameter 27ra/ A determines how many terms of this con¬ 
vergent series to be dominant [33] . 

In general, these scattering coefficients control the ex¬ 
trinsic states in the scattering field pattern. Once the 
composition of scatterer and geometry are given, the cor¬ 
responding scattering coefficients are determined. Con¬ 
versely, without knowing the composition in a scatterer, 
one can also have the same scattering coefficients. In the 
following, we introduce a systematic way to do the in¬ 
verse design for the scatterer by specifying the required 
scattering or absorption properties first. To do that, we 
introduce the phase diagram by expressing the scatter¬ 
ing coefficient as | exp{i6>i™’™^}, 

where the magnitude is a positive real value 

for every individual channel and is the corre¬ 

sponding phase shift. Since TM and TE modes are sepa¬ 
rable, we simplify the notations by neglecting the super¬ 
scripts, i.e., re-writing the amplitude square and phase 
of scattering coefficients as and 0^. In Fig. 1, by 

plotting the absorption cross-section of individual chan¬ 
nel, normalized to the common factor 27r/(2n + 1)A^, as 
a function of phase and amplitude square of scattering 
coefficient {0n^\Sn\‘^), we introduce the universal phase 
diagram for a passive scatterer. 

As shown in Fig. 1, even though we do not write down 
any exact formulas for this phase diagram, the range to 
support physical value for the amplitude square only ex¬ 
ists within 0 < ^ 1 ; while the phase is bound 

within 7r/2 < < 37r/2. It is worth to remark that 

this range to support our parameter space results from 
the intrinsic power conservation. We depict the allowed 
regions in colors. Outside the allowed regions, uncolored 
regions represent forbidden solutions for a passive scatter. 
Moreover, lossless absorption happens in the boundary 
between the allowed and forbidden regions. Along this 
lossless contour, there exists a family of solutions which 
all result in = 0, but can be supported with differ¬ 
ent scattering coefficients in amplitude and phase. It is 
known that with a localized surface plasmon in the sub¬ 
wavelength structure can assist lossless resonance con¬ 
dition EE], which corresponds to the point |S'^| = 1 



FIG. 1: (Color online) Phase diagram for each spherical har¬ 
monic channel is generated in terms of the phase On and am¬ 
plitude square for the corresponding normalized absorp¬ 
tion cross-section. Marked numbers shown in the contour 
lines correspond to the values of normalized absorption cross 
section of individual channels: 27rcr^^®/(2n-h 1)A^. Colored 
regions are physically allowed solutions; while uncolored re¬ 
gions (or in the White color) represent forbidden solutions. 
It is noted that the amplitude square is bounded within the 
range [0,1]; while the allowed phase is within [7r/2, 37r/2]. The 
Green cross-marker, localed at {On = 7T,\Sn\‘^ = 0.25), indi¬ 
cates the maximum value, 0.25, in the normalized absorption 
cross-section. 


and On = TT in our phase diagram. As for the invisible 
cloaks [HHIS], one can look for the supported solutions 
near the bottom of the phase diagram, i.e., \Sn\ =0, for 
dominant channels. 

Once the composited material in a scatterer has in¬ 
trinsic loss, the supported scattering coefficient moves to 
reside inside the colored region in this phase diagram. 
For each channel, the maximum value in the normalized 
absorption cross section is 27rcr^^®/(2n + 1)A^ = 1/4, i.e., 
the Green cross-marker shown in Fig. 1, corresponding 
to coherent perfect absorbers [MO], but which is also 
associated with the same amount of light scattered. The 
phase and amplitude of scattering coefficients to achieve 
a maximum absorption power is tt and 1/2, respectively. 
Moreover, along the contour for a constant absorption 
power, there exist a maximum and a minimum values in 
the scattering amplitude, both at the phase 0 = it. It 
implies that one may design a scatterer possessing the 
same absorption power, but with different scattering sig¬ 
nals. As for the concept on cloaking a sensor m to 
exact more information from the outside (correspond¬ 
ing to light absorption) but with the overall scattering 
performance being suppressed, in the phase diagram one 
can find the corresponding solutions located nt 0 = tt 







3 





FIG. 2: (Color online) Supported trajectories in the phase 
diagram are shown for different sets of the parameters: a, [3 
defined in Eq. (3). Here, trajectories with a constant 13 are 
shown in dotted-curves; while trajectories with a constant a 
are shown in dashed-curves. Two contours for a constant 
absorption power are also depicted in the Black color. 



FIG. 3: (Color online) Absorption and scattering cross- 
sections correspond to the contour shown in Fig. 2, which 
are depicted in terms of the parametric variable t defined in 
Eqs. (4-5). Here, the absorption power is a constant as re¬ 
quested, = 0.1; while there is a degree of freedom in the 
scattering power The insect illustrates the core-shell 

scatterer used as an example to design a passive electromag¬ 
netic devices with the constant absorption power. 


In particular, for the contour with a constant absorp¬ 
tion power in the phase diagram, as shown in the Black 
color in Fig. 2, one can use the parametric representation 
of the curve to describe this supported trajectory. The 
corresponding parameter sets of a and P can be found as 


(conjugate-matched condition) with the minimum scat¬ 
tering amplitude m- Even though above illustrations 
are demonstrated for a single channel, by including mul¬ 
tiple channels, the intrinsic single channel limitation can 
be broken to generate superscattering or superabsorber 
phenomena [30]. 

Through above examples, the phase diagram provides 
a universal map to all unusual scattering characteristics. 
Moreover, through material dispersion or the geometry in 
particle configurations, one can manipulate the allowed 
trajectories in the phase diagram to design a passive scat¬ 
terer with the required scattering and absorption prop¬ 
erties. For this purpose, we express the scattering coef¬ 
ficient for an isotropic and homogeneous scatterer with 
multiple layers as 

I iio'n iPn)^ ^ ^ 

^ n 

where Un and Vn are determinants of a 2N x 2N matrix 
constituted by spherical harmonic functions expressed at 
each boundary, and N is the number of layers in a single 
scatterer [3 nuns]. Here, we also rewrite this scatter¬ 
ing coefficient by introducing two real numbers: an and 
/3n, as shown in Eq. (3). By substituting Eq. (3) into 
the amplitude square and related phase of scattering co¬ 
efficients, one can have different trajectories in the phase 
diagram, as shown in Fig. 2. Supported trajectories for 
the parameter sets with a constant value of a or are 
plotted in dashed- and dotted-curves, respectively. 


an{t) 



(4) 


I3n{t) 



— cos{t), 

Qn 


(5) 


where the independent variable, t, is used for the para¬ 
metric equation, which is also bounded within the range 
t = [0,27r]. Here, for a given normalized absorption 
power, Qn = ( 2 n+T)A 2 ^n^^ oach channel, from Eqs. (4- 
5) one can see that the trajectories of an and Pn are 
following an elliptic equation. In terms of this paramet¬ 
ric variable t, in Fig. 3, we show the corresponding ab¬ 
sorption and scattering cross-sections for such a contour 
in the phase diagram. One can see that indeed the re¬ 
sulting absorption power is a constant as requested; 
while there is a degree of freedom to support different 
scattering powers 

In general, the parameter sets an and pn can be ma¬ 
nipulated through the geometry, material property, or 
the working frequency for a scatterer. As an example, 
we consider a passive scatterer in the configuration of a 
core-shell sphere, as illustrated in the insect of Fig. 3, 
which is composed by two concentric layers of isotropic 
and homogeneous materials. The geometrical parameters 
and material properties for this core-shell scatterer are 
the radius of core, Uc, the radius of whole particle, a, and 
^s{l^s)/ ^c{l^c) for the permittivity (permeability) in the 
shell/core regions, respectively. If the electrically small 
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FIG. 4: (Color online) The permittivities to support a con¬ 
stant absorption power are shown as a function of the para¬ 
metric variable t. For a given material in the shell region, 
Cs = 3.12, found solutions for the real and imaginary parts of 
the permittivity in the core region denoted as e'c and e'c are 
shown in (a) and (b), respectively. For a given material in 
the core region, Cc = 5, two families of found solutions shown 
in Eq. ( 8 ) are denoted as ef and e~ for the shell region. 
The corresponding real and imaginary parts of the permit¬ 
tivity denoted as e'g and e'/ are shown in (c, e) and (d, f), 
respectively. Results obtained from analytical formulas are 
depicted in solid-curves; while exact solutions from Mie’s the¬ 
ory are depicted in dashed-curves. In all cases, the core-shell 
geometries are all fixed with a = 1/24A and 7 = 0.9. 


approximation is satisfied for such a core-shell scatterer, 
it is known that the main contribution on the correspond¬ 


ing cross-sections dominantly comes from the dipole-wave 
scattering, i.e., n = 1. The requirement to satisfy is 
that all the size parameters are small enough, that is 
27ralX 1, 2'Ka-^esiJis/X 1, and 27Tay/e^jI^/X ^ 1 
for the outside, shell, and core regions, respectively. For 
such a two-layered scatterer, the corresponding scatter¬ 
ing coefficients are conducted from a 4 x 4 matrix by 
tracking TE and TM modes. However, for the electric 
dipole, the TM mode (n = 1) is dominant. By applying 
the continuity of electric magnetic fields established at 
the each boundary of shell-environment and core-shell, 
one can approximately express the coefficient /Ui^ 
as 

yTM ^ 2 A^ 27^(1 — es)(ec — Cg) — (2 + es)(ec + 2 es) 

U™ ~ 2 ( 2 ^a )3 73(e, - ee)(2e, + 1) + (1 - e,)(ee + 2e,) ’ 

(6) 

where 7 is defined as the ratio between the core radius 
to the whole particle radius, 7 = ac/a. If one replaces e 
by then we can obtain the other coefficient 
which for non-magnetic media is automatically zero for 
= Me = 1- By taking 7 = 1 or = Cc, result 

shown above can be reduced to the electric dipole equa¬ 

tion for a solid sphere. We want to mark that even though 
the electrically small approximation is applied here, how¬ 
ever, direct numerical results obtained by solving all the 
infinite channels only give a derivation around 10%. 

Now, for our core-shell system with the geometric size 
fixed, we provide a systematic way to find the correspond¬ 
ing material properties with a constant absorption power, 
as specified by the contour in the phase diagram shown 
in Fig. 2. To give a clear illustration, one may fix the 
material property in the shell or in the core region first. 
If we assume that the composition for the shell region 
is given, i.e., is fixed. Then, based on Eqs. (4-6), the 
corresponding solution for the permittivity in the core 
region is 


3(-2e, - 4 - 273 + 2^^es) - 2(ai + /?i)(27ra/A)3(2 - 2e, + 2e,73 + 7 ^) 

3(es + 2 + 2'y^€s - 2j^) + 2{ai + ^i)(27ra/A)3(l - e* - 26*73 - 


Solutions obtained from the analytical results shown in 
Eq. (7) are shown in Fig. 4(a) and 4(b) for the real and 
imaginary parts of the permittivity in the core region, re¬ 
spectively. In terms of the parametric variable. A, we can 
have a wide rang in selecting materials, and all of them 
have the same absorption power. A comparison between 
Fig. 3 and Figs. 4(a-b), we find that when A = 7r/2 the 
scattering power reaches a maximum value; while the 
for the required material has a minimum value. The rea¬ 
son comes from that dissipative loss is proportional to 
where E is the local electric field. In this sce¬ 
nario, with the help of a strong electric field, the loss to 


maintain the same absorption power can be minimized 
simultaneously. In Fig. 4, we also add the numerical 
calculations obtained by the exact solution from Mie’s 
theory, shown in the dotted-curves accordingly. A good 
agreement between our analytical solutions to the exactly 
numerical ones is demonstrated. 

On other hand, if the material property in the core 
region is specified, based on Eqs. (4-6) the corresponding 
solutions for the permittivity in the shell region to 
support a constant absorption power are governed by 

± _-9± ^/g'^ -4./h 
^ 1 


( 8 ) 
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with the short handed notations: 

/ = 2(l-73)[3-2(ai+/3i)(27ra/A)3], (9) 

g = 2 (ai + /3i)(27ra/A)3[73(l - 2e,) + 2 - e^] (10) 
+3 ( 27 ^ + 2'y^ec + + 4), 

h = ee(l-7")[6 + 2(ai+/3i)(27ra/A)3]. (11) 

There exists two families to support the materials in the 
shell region, denoted as e^. Again, we show the real and 
imaginary parts of the permittivity in the shell region 
for these two families in Fig. 4(c-d) and 4(e-f), as well 
as the exact solutions in dashed-curves, respectively. To 
one’s surprise, there exists a variety of choices for the 
material properties even some specific light scattering or 
absorption properties is given in the beginning. We want 
to emphasize that without the introduction of the phase 
diagram, it is not only hard to find the required mate¬ 
rial properties, but also difficult to recognize the power 
competition and limitation among these cross sections for 
each channel. 

Before the conclusion, we remark that when the electri¬ 
cally small approximation is not valid, multiple channels 
for TE and TM modes may be excited as expected. In 
this scenario, one can also apply our phase diagram, but 


not for a single channel only. Scattering coefficients from 
several dominant channels are just a natural extension by 
considering all of them onto the phase diagrams simulta¬ 
neously. Through the proper assignment to tune each 
dominant channels in the characteristic cross-sections, 
one may also have unique devices beyond the response 
of a single channel, such as subwavelength superscatter- 
ing and superabsorber pTH?!! |30]. 

In summary, in terms of the amplitude square and 
phase of the scattering coefficients, we introduce a univer¬ 
sal phase diagram as a compact tool to link the scatter¬ 
ing and absorption powers for each individual channel. 
Intrinsically, the power conservation gives the physical 
boundary in the parameter space for a passive scatterer. 
Not only the known exotic scattering and absorption phe¬ 
nomena can be illustrated in this diagram, supported tra¬ 
jectories are also demonstrated to design extrinsic-field- 
controllable scatterers. With the core-shell scatterers as 
an example, we reveal a systematic way to find a vari¬ 
ety of solutions in the composited materials, all of which 
possess the same absorption power. With the analogy 
among wave phenomena, the concept of our method can 
be ready applied to acoustic systems as well as quantum 
scattering system. 


[1] G. Mie, Annal. Phys. 330 377 (1908). 

[2] J. A. Schuller, E. S. Barnard, W.-S. Cai, Y. C. Jun, J. S. 
White and M. L. Brongersma Nat. Mater. 9, 193 (2010). 

[3] M. I. Tribelsky and B. S. Lukyanchuk, Phys. Rev. Lett. 
97, 263902 (2006). 

[4] F. Monticone, C. Argyropoulos, and A. Alu, Phys. Rev. 
Lett. 110, 113901 (2013). 

[5] A. Alu and N. Engheta, J. Appl. Phys. 97, 094310 (2005). 

[6] R. E. Hamam, A. Karalis, J. D. Joannopoulos, and M. 
Soljacic, Phys. Rev. A 75, 053801 (2007). 

[7] H. Noh, Y. Chong, A. D. Stone, and H. Cao, Phys. Rev. 
Lett. 108, 186805 (2012). 

[8] H. Noh, S. M. Popoff, and H. Cao, Opt. Express 21, 
17435 (2013). 

[9] M. 1. Tribelsky, Europhys. Lett. 94, 14004 (2011). 

[10] V. Grigoriev, N. Bonod, J. Wenger, and B. Stout, ACS 
Photonics 2, 263 (2015). 

[11] A. Alu and N. Engheta, Phys. Rev. E 72, 016623 (2005). 

[12] A. Alu and N. Engheta, Opt. Express 15, 3318 (2007). 

[13] S. Muhlig, M. Farhat, C. Rockstuhl, and F. Lederer, 
Phys. Rev. B 83, 195116 (2011). 

[14] S. Muhlig, A. Cunningham, J. Dintinger, M. Farhat, S. B. 
Hasan, T. Scharf, T. Burgi, F. Lederer and C. Rockstuhl, 
Sci. Rep. 3, 2328 (2013). 

[15] M. Farhat, S. Muhlig, C. Rockstuhl, and F. Lederer, Opt. 
Express 20, 13896 (2012). 

[16] A. Alu and N. Engheta, Phys. Rev. Lett. 102, 233901 
(2009). 

[17] Z. Ruan and S. Fan, Phys. Rev. Lett. 105, 013901 (2010). 

[18] Z. Ruan and S. Fan, Appl. Phys. Lett. 98, 043101 (2011). 

[19] A. Mirzaei, A. E. Miroshnichenko, 1. V. Shadrivov, and 
Y. S. Kivshar, Appl. Phys. Lett. 105, 011109 (2014). 

[20] N. M. Estakhri and A. Alu, Phys. Rev. B. 89, 121416(R) 


(2014). 

[21] Z. Ruan and S. Fan, J. Phys. Chem. C 114, 7324 (2010). 

[22] K.R. Catchpole and A. Polman, Opt. Express 16, 21793 
(2008). 

[23] X. Zhang, Y. L. Chen, R.-S. Liu, and D. P. Tsai, Rep. 
Prog. Phys. 76, 046401 (2013). 

[24] M.-C. Tsai, J.-Y. Lee, P.-C. Chen, Y.-W. Chang, Y.-C. 
Chang, M.-H. Yang, H.-T. Chiu, I.-N. Lin, R.-K. Lee, 
C.-Y. Lee, Appl. Catal. B 147, 499 (2014). 

[25] L. R. Hirsch, R. J. Stafford, J. A. Bankson, S. R. Sershen, 
B. Rivera, R. E. Price, J. D. Hazle, N. J. Halas, and J. L. 
West, Proc. Natl. Acad. Sci. (PNAS) U.S.A. 100, 13549 
(2003). 

[26] M. 1. Tribelsky, A. E. Miroshnichenko, Yu. S. Kivshar, 
B. S. Lukyanchuk and A. R. Khokhlov, Phys. Rev. X 1, 
021024 (2011). 

[27] J. W. Shi, F. Monticone, S. Elias, Y. W. Wu, D. Ratch- 
ford, X. Q. Li, and A. Alu, Nat. Commun. 5, 3896 (2014). 

[28] Y. Pu, R. Grange, c-L Hsieh, and D. Psaltis, Phys. Rev. 
Lett. 104, 207402 (2010). 

[29] S.-W. Chu, T.-Y Su, R. Oketani, Y.-T. Huang, H.-Y. 
Wu, Y. Yonemaru, M. Yamanaka, H. Lee, G.-Y. Zhuo, 
M.-Y. Lee, S. Kawata, and K. Fujita, Phys. Rev. Lett. 
112, 017402 (2014). 

[30] A. Mirzaei, 1. V. Shadrivov, A. E. Miroshnichenko, and 
Yu. S. Kivshar, Opt. Express 21, 10454 (2013). 

[31] R. Fleury, J. Soric, and A. Alu, Phys. Rev. B. 89, 045122 
(2014). 

[32] M. Born and E. Wolf, “Principle of optics,” 7th Ed. 
(Cambridge University Press, 1999). 

[33] C. F. Bohren and D. R. Huffman, “Absorption and scat¬ 
tering of light by small particles,” (Wiley, 1983). 

[34] S. A. Maier, “Plasmonics: fundamentals and applica- 



6 


tions,” ( Springer, 2007). 



